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Abstract 

We study the scattering of atoms, rotons and phonons at the free surface of 
4 He at normal incidence and calculate the evaporation, condensation and re- 
flection probabilities. Assuming elastic one-to-one processes and using general 
properties of the scattering matrix, such as unitarity and time reversal, we 
argue that all nonzero probabilities can be written in terms of a single energy- 
dependent parameter. Quantitative predictions are obtained using linearized 

time dependent density functional theory. 
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Superfluid helium at low temperature exhibits the peculiar phenomenon of quantum 
evaporation. Elementary excitations, such as rotons and high energy phonons, propagate in 
the liquid with long mean free path and with an energy which can exceed the binding energy 
of atoms in the liquid. Thus, when an elementary excitation impinges upon the free surface it 
may eject an atom out of the liquid through a quantum process. The scattering of excitations 
and atoms at the free surface of the superfluid have been extensively studied in the recent 
years (see for instance Refs. HHI! and references therein), but the comparison between theory 
and experiments is still unsatisfactory. Previous theoretical approaches to the problem of 
quantum evaporation made use of semiclassical approximations 0J|. Recently |||7j we 
have calculated the probability of evaporation, condensation and reflection of rotons and 
atoms using a time dependent density functional (TDDF) theory. The density functional 
has proven a powerful method to investigate various structural features of inhomogeneous 
superfluid helium || including static and dynamic properties of the free surface, droplets and 
films. When applied to quantum evaporation, the linearized TDDF theory accounts only for 
one-to-one processes, but includes quantum effects beyond the semiclassical approximation. 
So far we have restricted the analysis to the phonon forbidden region |||7|], corresponding 
to incidence angles such as phonons are excluded by energy and momentum conservation. 
In the present work we apply the same theory to the case of scattering at normal incidence, 
where atoms, rotons as well as phonons take part in the scattering processes. We show that 
the complexity of the scattering processes is dramatically reduced on the basis of simple 
arguments based on symmetry properties of the scattering matrix and threshold effects. 
The final result is that all nonzero probabilities can be expressed in terms of a single energy- 
dependent parameter. The numerical solution of the time dependent density functional 
equations confirms this picture and provides quantitative predictions for the probabilities as 
a function of energy. 

In Fig. 1 we show the phonon-roton dispersion curve in bulk liquid. The phonon branch 
goes from q = up to the maxon (A*). On the right of the maxon the slope is negative 
and corresponds to the dispersion of R~ rotons, that is, rotons with negative group velocity. 
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Then the slope becomes positive again for R + rotons on the right of the roton minimum 
(A). The threshold for atom evaporation at normal incidence coincides with the chemical 
potential, \fx\ = 7.15 K (dashed line): only excitations with energy higher than can 
produce quantum evaporation. The energy of an evaporated atom is simply + h 2 q 2 /2m. 
For scattering at normal incidence, the problem is unidimensional; the wave vector q is 
orthogonal to the free surface and is not a conserved quantity, since the surface breaks 
translational invariance. As already said, we consider elastic processes in which the energy 
Hcu is conserved. 

Four types of excitations interact at the free surface of the superfluid: atoms (a), phonons 
(p), rotons with positive (+) and negative (— ) group velocity. One can define the scattering 
matrix by means of ^ out = S^ in , where ^/ in and ^ out are the incoming and outgoing 
asymptotic solutions, respectively. The matrix element connects a particular input (i) 
and output (j) channel, with i,j = a,p, +, — , and depends on energy and incident angle. 
It determines the probability associated with each scattering process by = \Sij\ 2 . These 
channels correspond to atoms, phonons and rotons excitations states. Thus, there are 16 
complex scattering matrix elements S^, and hence 16 probabilities P^-, to be determined. 
The unitarity and the time reversal symmetry (t — > —t) of S imply = S^ 1 and S* = S' 1 , 
respectively. Combining these two conditions one finds that the scattering matrix elements 
must satisfy the general property SV,- = Sjt, reducing the independent matrix elements to 10. 
Furthermore, the unitarity condition S*S = 1 gives ten additional constraints. For instance, 
if the incident excitation is of type i, then one has the unitarity conditions 

£Kf = i, (i) 

where i, j = a,p, +, — , yielding four linear combinations of matrix elements. The other six 
conditions are 

Saa^ap + S* p S pp + S*_S- p + S* + S +p = (2) 

SLSa- + s* ap s p _ + s:_s^ + s: + s + ^ = (3) 
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Saa,Sa+ + S* p S p+ + S*_S- + + S* + S ++ — (4) 

S pa S a - + S pp Sp^ + S p _S — + S p+ S^ — = (5) 
Spa,Sa+ + S pp Sp+ + S p _S- + + S' p+ S' ++ = (6) 

S_ a S a + + S_ p S p + + S 1 5 — |- + S = . (7) 

We will now combine these rigorous relations with further arguments about threshold effects 
and momentum exchange at the surface in order to reduce the number of relevant unknowns. 

The roton minimum represents the threshold energy to excite R~ and R + rotons in the 
liquid. These excitations coincide at A and, on the basis of symmetry arguments, one can 
prove that the mode-change reflection between R~ and R + is the dominant one just above 
A (P+- — 1) and that the other probabilities involving R~ and R + are small. Conversely, 
below it only phonons and atoms are present. Among the different processes the phonon <-> 
atom scattering is expected to be favoured (P pa — 1), since it implies the smallest change 
of momentum. In fact, the normal reflection of a phonon or an atom would imply a much 
larger momentum transferred to the surface and this seems unlikely for a smooth surface like 
the one of helium, except near the threshold hu ~ where the atom momentum goes to 
zero and full reflection takes place. Now, assuming P + _ ~ 1 and P pa ~ 1 into the unitarity 
conditions (d), one finds that all the remaining probabilities should vanish close to A. 

Similarly, the maxon is a threshold for phonons and R~ rotons and one expects the 
mode-change reflection between R~ and phonons to dominate just below A* (P p - — 1). 
Just above it, one has only R + and atoms. For the same argument of smallest momentum 
transferred to the surface, the R + <->• atom scattering should be much favoured (P+ a — 1) 
with respect to the normal reflection of rotons and atoms. Putting P p _ ~ 1 and P +a ~ 1 in 
the unitarity conditions one finds again that the other probabilities have to vanish near A*. 

The above arguments suggest that P pa and P + _ decrease from 1 to by increasing the 
energy from A to A*, while P p _ and P +a increase form to 1 in the same range. All 
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probabilities should be smooth functions of the energy. The ones that are both at A and 
A* are expected to be small everywhere in between. If we make the assumption that they 
vanish for any value of energy (A < Hlu < A*), 

Paa — Ppp — P — P++ — P-a = P+p = , (8) 

we find simple and useful relations among the remaining nonzero probabilities. In fact, by 
inserting assumption (||) in the unitarity conditions (IU-0), after some simple algebra one gets 


P pa = P + _ = 1-P p _ = l- P +a . (9) 

This means that all nonzero probabilities at normal incidence can be written in terms of 
a single parameter. This result is expected to hold for any theory accounting only for 
one-to-one elastic processes. 

In the remaining part of the paper we show that the numerical calculation of Pij, within 
time dependent density functional theory, are consistent with the general properties of the 
scattering matrix and with the threshold effect near A and A*. Furthermore, it confirms 
the validity of assumption (g). 

In order to calculate explicitly the probabilities one needs a theory for the elementary 
excitations of the nonuniform liquid. The theory must provide a spectrum of elementary 
excitations close to the experimental one, since the proper energy balance between the exci- 
tations is crucial; moreover, it must allow one to calculate the asymptotic flux of elementary 
excitations in a given process, in order to extract the corresponding matrix elements. 

In a density functional approach |],[7j one assumes the energy of the system to be in 
the form E = J dr7i[^/, where the complex function \I/ = <& exp(iS/h) is related to the 
density and velocity of the superfluid by means of p = $ 2 and v = (l/m)VS. We use 
a phenomenological energy density 7i which has been adjusted to provide an accurate 
description of the equation of state, the static response function and the phonon-roton 
dispersion law of the bulk liquid. The description of the functional can be found in our 
previous papers EH- 
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The eigenenergies and eigenf unctions of the system are calculated within a time depen- 
dent density functional approach. The propagation of the elementary excitations is described 
by linearizing the wave function around the ground state \l/o( r ) : 

#(r,*) = ¥ (r) + <Stt(r,t). (10) 

Far away from the surface S^(r,t) corresponds to the propagation of phonons and rotons 
in the liquid and free atoms in the vacuum. Due to linearization, only one-to-one processes 
are taken into account. Thus the theory can not describe inelastic processes, such as decay 
into multi-ripplons or multi-phonons. 

For practical reasons we work in a slab geometry. The semi-infinite system is simulated 
by slabs of liquid with thickness ranging from 50 to 100 A, and centred in a finite box 
whose size is of the order of 100 -j- 150 A. The density functional provides also the ground 
state density profile of the liquid as the stationary solution of lowest energy. The profile 
of the liquid-vacuum interface has a thickness of about 6 A and a shape close to the one 
obtained with ab initio Monte Carlo calculations (see Ref. || for details). Choosing z along 
the normal to the surface, one can write the wave function in the form 

5^(z,t) = f(z)e'^ + g(z)e^. (11) 

The quantities f(z) and g(z) are real functions which have to be determined, together with 
u), by solving self-consistently the equations of motion 

5 Jdtjdr **] - /zW - **^^| = , (12) 

linearized with respect to / and g. Equations (|TTD and ( 12]) assume the typical form of 
the random phase approximation (RPA) for bosons. In the formalism of RPA, f(z) and 
g(z) take into account the particle-hole and hole-particle transitions, respectively. The two 
components / and g of the excited states are coupled by the equations of motion (p~2j). 

We have solved numerically the linearized TDDF equations in the finite box, obtaining 
a set of discrete eigenenergies uj and the corresponding eigenf unctions f(z) and g(z). For 
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a given to, the solution is a stationary state. Inside the slab, f(z) and g(z) are oscillating 
functions associated with the propagation of phonons and rotons, whose dispersion law 
is shown in Fig. 1. Outside the slab, where particles are uncorrelated, the function g(z) 
vanishes while the equation for f(z) coincides with the Schrodinger equation for free atoms. 
By looking at the Fourier transforms of the signal far away from the surface, both inside and 
outside the slab, one can evaluate the asymptotic amplitudes (fi and gi) and the momentum 
qi of each type of excitation which contributes to 5^. A fitting procedure has been used in 
order to extract the numerical values of fi and gi needed for the analysis of the evaporation 
rates. 

A given scattering process at a certain energy can be obtained as a linear combination of 
different stationary solutions at that energy. The latter can be found by slightly varying the 
slab thickness (L s i a b) and the box size (Lf, ox ). Then, one can evaluate the flux of incoming 
and outgoing excitations. In the present linearized TDDF theory, the current associated 
with a given elementary excitation is 

ji = v<(|/,| 2 -|ft| 2 ), (13) 

where v is its group velocity. From the asymptotic fluxes one can easily evaluate the evap- 
oration, condensation and reflection probabilities 0. 

We have first verified that the resulting probabilities satisfy, within the accuracy of the 
calculation, the symmetry property = Pji and the unitarity conditions (|l|) in the whole 
range of energy hu > Explicit results are shown in Table |for an intermediate value of 
energy (%u = 11 K). Four states at the same energy are combined to describe the scattering 
processes (atom — > atom, phonon, R~ , R + ), (phonon — > atom, phonon, R~,R + ), (R + — > 
atom, phonon, R~ ,R + ) and (R~ — ► atom, phonon, R~,R + ). The numerical uncertainty, 
which is less than ±0.05 for all probabilities, originate mainly from the fact that the states 
in the linear combinations may be not enough linearly independent [[!]]. 

We find also that the probabilities P aa , P pp , P , P ++ , P- a and P p+ are extremely small at 

all energies; they turn out to be zero within the present accuracy, confirming the hypothesis 
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(§) made before. Furthermore, the nonzero probabilities turn out to verify relations @. 
It is therefore possible to resume all the information about the nonzero probabilities in a 
single function of energy, P{oj). Let us call P{oj) = P+ a - This quantity is plotted in Fig. 2, 
in the energy range A < Tuo < A*. The probability of mode-change scattering between 
R~ and phonons (P p -) is equal to the evaporation probability for R + rotons (P+ ); both 
probabilities start from at the roton minimum and increase to 1 at the maxon energy. 
The evaporation probability for R + rotons is equal to 1 even above the maxon energy. The 
evaporation probability for phonons (P pa ) is equal to the one for the mode-change process 
between rotons (-P+-). They are equal to (1 — P +a ), so that they decrease from 1 to between 
A and A*. It is worth noting that the present results for the evaporation probability P +a 
at normal incidence coincide, within the error bar, with the ones obtained previously in the 
phonon forbidden region N7J, that is, for R + rotons impinging at about 15 -r 20 degrees. The 



same is true for the mode-change probability P H 1 1 I 



Only few experimental estimates are available for the probabilities Py. Indeed, by mea- 
suring time-of-flight and angular distribution of evaporated atoms, one-to-one evaporation 
processes have been clearly seen [|IJ, but a quantitative determination of the ratio between 
incoming and outgoing fluxes, and hence P^-, is difficult. There are evidences for a sizable 
probability P+ a , which should increase with a trend similar to the one in Fig. 1 JTJj. A recent 
estimate ]13| of the phonon — > atom probability, for high energy phonons (Tuo > 10K), is 
P pa ~ 0.1. Our theory, in the same range of energy, gives a value which decreases rapidly 
from 0.5 to 0. We also find that an incident atom condenses with probability P ap + P a+ ~ 1 
and hence the reflection probability P aa is almost zero. This is in agreement with the ob- 
served small reflectivity but the experiments support also the idea that the atom 
condensation might occur through processes like two-phonon or multi ripplon production. 
The extension of the present formalism to include such mechanisms, beyond the one-to-one 
hypothesis, remains an important task, for a more systematic comparison between theory 
and experiments. 
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FIGURES 

FIG. 1. Phonon-roton dispersion in liquid 4 He. The dashed line is the threshold = 7.15 K 
for atom evaporation. 

FIG. 2. Probability P = P+ a as a function of energy. The other nonzero probabilities can be 
extracted using Eq (||). Energy scale starts from A. 
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TABLES 

TABLE I. Probabilities for four different scattering processes (incident atom, phonon, R~ 
and R + , respectively) described by linear combinations of four states at the same energy (11 K) 
and for different values of (L s i a ij,Lb ox ). The unitarity condition, Eq. (Q), is labeled by E. 



incident exc. 


P 


E 


atom 


P aa =0.0002 


P ap =0.3384 


Pa- 


.=0.0003 


Pa+ 


=0.6837 


1.0226 


phonon 


P pa =0.3236 


P pp =0.0020 




_=0.6487 


Pp+ 


=0.00002 


0.9743 


R- 


P_ a =0.0004 


P_ p =0.6859 


P__ 


-=0.0013 




.=0.3408 


1.0284 


R+ 


P +a =0.6537 


P +p =0.0001 


P+ 


_=0.3222 


P++ 


=0.00002 


0.9760 
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